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We study generalized Korteweg-de Vries (KdV) equations derivable from the Lagrangian: 
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The usual field u{x,t) of the generalized KdV equation is defined by u{x,t) = ipx{x,t). 
The equation of motion derived from this Lagrangian has solitary wave solutions of both 
the usual (non-compact) and compact variety ( "compactons" ) . For the particular case that 
p = m = n + I, the solitary wave solutions have compact support and the feature that 
their width is independent of the amplitude. We discuss the Hamiltonian structure of these 
theories and find that mass, momentum, and energy are conserved. We find in general that 
these are not completely integrable systems. Numerical simulations show that an arbitrary 
compact initial wave packet whose width is wider than that of a compacton breaks up into 
several compactons all having the same width. Upon scattering the compactons found here 
behave similarly to those found in the equations of Rosenau and Hyman, the scattering 
is almost elastic, with the left over wake eventually turning into compacton-anticompacton 
pairs. Often there are two different compacton solutions for a single set of parameters having 
the same generic form Acos^{d(,), where ^ = x — ci and —tt/2 < < tt/2. When this is the 
case, the wider solution is stable, and this solution is a minimum of the Hamiltonian. 
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I. INTRODUCTION 

The observed stationary and dynamical patterns in na- 
ture are usually finite in extent. However most equations 
that admit solitary wave solutions yield solutions that are 
infinite in extent, although of a localized nature. There- 
fore, the recently discovered solitary waves with compact 
support (compactons) found by Rosenau and Hyman 0] 
represents a welcome development. Rosenau and Hyman 
found in numerical experiments that the compact solitary 
waves on collision reemerge as compactons, with a tiny 
amount of energy going into a zero mass ripple which 
eventually reemerges into a compacton- anticompacton 
pair of small (< 5%) amplitude. The original equation 
studied by Rosenau and Hyman was ||l| 



Unlike the ordinary KdV equation the above equation 
is not derivable from a first order Lagrangian except for 
n = 1, and did not possess the usual conservation laws 
of energy, momentum, and mass that the KdV equation 
possessed. Because of this, we considered earlier a differ- 
ent generalization of the KdV equation based on a first 
order Lagrangian formulation. That is, we considered j^] 
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For 1 < m = n < 3 this equation had solitary wave 
solutions of the form [cos(a^)]^/'^™~^\ where ^ — x — ct 
and — 7r/2 < < 7r/2 . For m==2,3 they obtained: 



This Lagrangian leads to a generalized sequence of 
KdV equations of the form: 

K*{l,p) : Ut + u'-^^Ux + a{2uPuxxx + 'Ipu^^'^UxUxx 

^p{p-l)u^-\uxf) = Q. (1.4) 



where 



if(2,2):w, = |cos2(e/4) 



u{x,t) = Lpx{x,t). 
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These equations have the same terms as the equations 
considered by Rosenau and Hyman, but the relative 
weights of the terms are different. In that study we 
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showed that for < p < 2 and I = p + 2, these mod- 
els admit compacton solutions for which the width is in- 
dependent of the amplitude. Since higher order KdV 
equations also admit soliton solutions, it is interesting to 
enquire whether higher order generalized KdV equations 
also admit compacton solutions, and whether there are 
cases where the width is independent of the amplitude. 
This is the first motivation for the present study. The 
second motivation is to study the scattering properties 
of the compactons numerically to see if they behave sim- 
ilarly to those found by Rosenau and Hyman. The gen- 
eralization which we will study in this paper is described 
by the Lagrangian: 



(1.6) 



which includes one extra term with higher derivatives. 
This Lagrangian Eq. (1.6) generalizes both the usual 
KdV Lagrangian and our previous generalized KdV La- 
grangian, and shares with both equations that because 
of the invariance of the action under time, and space 
translations as well as the shift of the field by a constant 
{ip ^ ip + a) that energy, momentum and mass are con- 
served by Noether's theorem. 

The rest of the paper is organized follows: In section 2 
we will discuss the Hamiltonian structure of these classes 
of theories and determine the equation for the travelling 
wave solutions to (1.6). In section 3 we give arguments 
based on a variational approach (see |^ Q) as well as 
dimensional reasons for there being solutions with the 
width of the solitary wave being independent of the am- 
plitude and we determine the connection between the 
energy and momentum of these solitary waves: namely 
E oc Pc. In particular we find that the condition for hav- 
ing solitary waves with width independent of amplitude 
is p = m = n + l. In section 4 we obtain exact compacton 
solutions of the form A cos^ d{x — ct) and verify that the 
relations among the global variables obtained using the 
variational approach are exact. In section 5 we give nu- 
merical results for the scattering of two compactons, as 
well as for the breakup of an arbitrary wave packet into 
several compactons. In an appendix we show that most 
of the exact solutions as well as the criteria for their sta- 
bility could be obtained by assuming the exact compact 
ansatz and minimizing the action on that class of func- 
tions. The unstable solutions are maxima of the effective 
Hamiltonian. 



II. 



GENERALIZED KDV EQUATION AND 
PROPERTIES 



From the Lagrangian Eq.( 1.6) we obtain the general- 
ized KdV equation: 



a 



p+l 
7n 



,m-l„,2\ 



+ jd^xxiu^'ulu^x) = 0, 
as well as the conserved Hamiltonian 
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where u{x,t) — (p xjx, t). We notice that the Lagrangian 
given by Eq. ( [L.6| ) is invariant under the transformations 
(i) f{x,t) ip(x^t)+ci; {u)x X+C2 and (iii) t t+ca, 
where ci, C2 and C3 are constants. By a direct application 
of Noether's theorem this leads to the three conservation 
laws of mass M, momentum P and energy E, where E 
is given by Eq. (2.2) and M and P are given by 
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u^{x,t)dx. (2.3) 



The equation of motion Eq. ( |2.l[ ) is also invariant under 
the transformations as given by 



provided that 
n + I 
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However, as seen below for the interesting compacton 
case where m = p = n + I, these equations are invariant 
under the more restricted transformations: 



ku^ X ^ x^ ij ' ^ k ^ ij . 
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We shall see later that this condition is necessary so as 
to have compactons whose width is independent of the 
amplitude. As an aside we would like to point out that all 
the compactons so far discovered with width independent 
of amplitude have the p rope rty that the field equations 
are invariant under Eq.( 2.4). 

The canonical structure of these theories is similar to 
that found in in that by postulating that the u{x) 
satisfy the Poisson bracket structure H 



{u{x),u{y)} = dx5{x - y). 
we obtain that 

^ SH ^ 



(2.5) 
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with H being given by Eq. ( . We also find that with 
our definition of P given by Eq.( pT^ ), that P is indeed 
the generator of the space translations: 



(2.7) 
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Since our equation is a generalization of the equation dis- 
cussed in y], one is also not able to show the existence 
of a bi-Hamiltonian structure using the conserved mo- 
mentum as a possible second Hamiltonian as was done 
for the ordinary KdV equation. So on these grounds one 
expects that our general Lagrangian does not correspond 
to an exactly integrable system, except for the original 
KdV equation case. This will be verified by our numeri- 
cal results on scattering where there is some energy going 
into compacton pair production following scattering. 



A. Equation for solitary waves 



If we assume a solution to in the form of a trav- 
elling wave: 



u{x,t) = f{x- ct) = f{y), 
one obtains for f{y): 

rrn-l r2\ 
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Integrating once we obtain: 
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where ci is a constant of integration. This equation needs 
several more integrations before a solution in terms of 
quadrature is obtained, unlike the previous equation we 
studied where two integrations were sufficient. Thus an 
explicit solution in terms of quadratures is not available 
and one must use "educated" guess ansatze to find exact 
solutions. It is therefore quite useful to look at a simple 
variational approach to the solitary wave problem to un- 
derstand some of the global features of the solitary waves. 
Previous experience |Q] |^] with variational methods have 
shown that many of the global relations among conserved 
quantities as well as conditions for the width to be inde- 
pendent of the amplitude can be obtained from a simple 
variational ansatz. Also whether the solitary wave is ex- 
pected to be stable or not can be inferred from the fact 
that the variational solution is a minimum of the action 
or just a saddle point of the action. Before addressing 
these questions, we just remind the reader that at 7 = 
the Lagrangian we are studying reduces to the previous 
generalized KdV problem we studied |^ so that these 
equations include all the KdV solitons and compactons 
that we discussed earli er as a special case. We will find 
exact solutions of Eq. ( 2.1C ) by inserting the ansatz: 



f{y) = Acqs'- (dy) 



(2.11) 



and determining the parameters A, r, and d by a consis- 
tency argument. 



III. VARIATIONAL APPROACH 

Our time-dependent variational approach for studying 
solitary waves is based on the principle of least Action. In 
previous work Q [|j 0, we introduced a post-Gaussian 
variational approximation, a continuous family of trial 
variational functions more general than Gaussians, which 
can still be treated analytically. We assumed a varia- 
tional ansatz of the form 
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i„(x,t) = yl(t)exp -b{t)\x - q{t)\ 



The variational parameters have a simple interpretation 
in terms of expectation values with respect to the "prob- 
ability" P 
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where the conserved momentum P is defined as above 
P=^J[u,ix,t)fdx. (3.2) 

We have that 
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Extremizing the effective action for the trial wave func- 
tion Uy leads to Lagrange's equations for the variational 
parameters. We find that for all values of the parame- 
ters (1,171, n,p) the dynamics of the variational parame- 
ters lead to solitary waves moving with constant velocity 
and constant width. However the solutions found are of- 
ten maxima or saddle points of the action. We will find 
that when this happens the exact solitary wave solution 
is unstable. For the special case of p = m = 71 -I- ^ we 
will show below that the width of the soliton is indepen- 
dent of the amplitude and velocity. For that case we also 
obtain the relationship: 



E = -- 
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We will find that the exact solitary wave solutions satisfy 
this relati onship as long as the integration constant ci of 
Eq. ( |2.10D is zero. 

The starting point for the variational calculation is the 
action 



Ldt, (3.5) 



where L is given by ( |1.3| ). 

Inserting the trial wave function u„ we obtain: 

T{q, (3, P,s)= J dt[-Pq - Heff]; (3.6) 

where H^f / is the Hamiltonian evaluated using the vari- 
ational wave function Uy. We find 

H,ff = Ci(s)6^'/4^p(f+2)/2 + C2(s)6(™+4)/4-f("+2)/2 

+ C'3(s)&(«+3'+8)/4sp(«+i+2)/2_ ^3 

where Ci, C2 and C3 are functions of s, a, (3, 7 but inde- 
pendent of P, b. We notice that there is no momentum 
conjugate to b. We eliminate the variable of constraint b 
(using 5T/db = 0) and obtain the equation: 



0= ^Ci(s)dfp(f+2)/2 + !^ 
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From Eq. (3 



d = b^. 
we see that when 

p = m ^ n + I 



(3.8) 



(3.9) 
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the momentum P factors out of the equation, so that 
the width of the soliton which depends on b does not 
depend on P and thus is independent of the amplitude or 
velocity. This special case is precisely the case when the 
exact solution is a compacton whose width is independent 
of amplitude as we shall show below. From Lagrange's 
equation 5T/5P = we have that 
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When p = m = n + I this yields the relationship: 
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(3.12) 



which we will verify is true for the exact compacton solu- 
tions when ci = 0. For this special case p — m = n + l 
the functions Ci are given by: 
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[p + l){p + 2)s{p + 2)^ 



Cs = ~{2sy+\-iyN^{p + 2) 

(2s - l)[p^(2s - 1) -I- 2p(l + I - 2ls) + f{2s - 1) 
+ 2/(s-l) + 4s-l] (3.13) 

where iV = 2^/[r(l-K^)]'^. 

For the special case that p — m = n + I our procedure 
for determining the optimal trial wave function is very 
straight forward. The effective Hamiltonian simplifies to 
(from now on we drop the subscript on Hgff): 

Hid,s) = P^dP[Ci{s)+d*C2{s) + d^'+^C3is)] (3.14) 

Minimizing with respect to the width parameter d, we 
obtain: 

pCi + {p + 4)d*C2 + {p + 2l + 8)d2'+8C3 = (3.15) 

For I ~ 0, one can analytically solve this for d{s) and one 
usually finds that there are 2 real solution in the vicin- 
ity of s = 1 which is the Gaussian trial wave function. 
Inserting these 2 solutions {di{s)) into H we obtain an 
analytic expression of the form 

H,is) = H{d,{s),s) 

Imposing dHi/ds — 0, we find that the two solutions 
{di(si), Si\ (i=l,2) have the property that one is a saddle 
of H in the {d, s} plane and one is a minimum. We find 
in our numerical simulation that only the solutions that 
is a minimum is stable. In general it is easy to find the 
extrema of H in the {d, s} plane graphically by making 
a 3 dimensional plot of H . What is convenient about the 
class of trial wave functions discussed here, is that one 
can explicitly determine the effective action for all values 
of p, I, m, n and get the general features of the results for 
all values of the parameters. In the appendix we will use 
the form of the exact compacton wave function in the 
variational approach. By doing so we will recover most 
of the exact compacton solutions that we find in the next 
section. We also obtain the result that when there are 
multiple solutions the stable one is a minimum of the 
effective Hamiltonian for the exact shape function. 



IV. EXACT COMPACTON SOLUTIONS 

Motivated by our variational results we will concen- 
trate on the particular case where p^m = n + l^Q. 
which corresponds to Compactons whose width is inde- 
pendent of amplitude. Assuming a solution of the form 



— 7r 7r 
u{x,t) — Acos^[d{x — ct)]; < d{x — ct) < — 

TT 

u{x,t) = 0; \d{x — ct)\>—, 



(4.1) 
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we look for consistent solutions for A, r,c and d in terms 
of a, P and 7 and ci which are real. We check whether 
there is an analogous variational solution of the post- 
Gaussian variety corresponding to a minimum and not 
just a saddle point or maximum of the action. We also 
show in the appendix that if we a assume a trial wave 
function of the exact type: 

u{x,t) = Acqs'' [d{t){x - q{t)] 

the exact solutions are either maxima or minima of the 
effective Hamiltonian as a function of d. We find that 
when there is a unique solution, then it is a minimum 
of the Energy. When there are two solution, the wider 
one is stable and is also a minimum of the Energy. The 
narrower one is unstable and is a maximum of the Energy. 
Now let us look at some special cases. 



A. p - 



n and / = case 



For the case p = m ~ n and Z = it is possible to find 
a general class of solutions for arbitrary p. Inserting a 
trial solution of the form: 

u{y) = A cos'' {dy) (4.2) 

into Equation( ^.lO ) we obtain the consistency equation: 



0= -3ci + 2Acx'' 
+A^+Px-^+''+''P (1 - r) rd'^j X 
(12 - lOr + 2r^ - llrp + 5r^p + 2rV) 



-2^1- 



„-2+r 



+''Prd^ (-2 + 2r + rp) x 



(-/? + 2^27 - 2rd'^'^ + r^d^-1 - rd^-fp + 2r'^d^^p) 
+A^+Pr'^d^ (2 + p) (2/3 - r'^d^-1 - 2r^d^-fp) x''+''p 
2 



l+p 



■A^+Pax 



r+r p 



Here x = cos{dy). All the powers of x must have zero 
coefficient for the trial solution to be an actual solution. 
This leads to various conditions depending on the values 
of r and p. If r{l+p) = 4 then there can also be solutions 
with Ci being nonzero. First let us consider the case when 
ci = 0. In that case for consistency we need either rp = 2 
or rp = 4. 



1. rp=2 



The first condition tells us that either 7 = or r = 1. 

When 7 = 0, we get the solution we found in our earlier 
work 0: Namely: 



ap 



4/3 a+p) i2+p) 



^ c(l+p) (2+p) 
2 a 

When 7 7^ we instead get the solution r ~ 1 {p — 2) 
and 

A^^ 



2 /3 d2 _ 6 d4 7 
and two possible solutions for the width: 



^2 _ 12/3± ^144/32 - 120^7 
~ 6O7 

Which means we also can write the equation for A"^ as 

c= [a-2(3d^) 

o 

Thus when 7 7^ one only gets a solution when p = 2. 
A particular case of this solution is: a = 6, 7 = 3 and 
/3 = 4. Then there are two solutions having — 1/3 or 
^2 ^ 1/5, The first is 



u = ^3c/2cos( ^ 

v3 

The conserved quantities for this solution are 



(4.5) 



M = 3\/2Z; P = -\/3c7r; E 



16 



(4.3) Thus we find that the relationship: 



E 
P 



2c 
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that we derived from our variational approach is exact 
here. 

The second solution : 



/25c y 

u = \ cos — — 

V 22 V5 

has as its conserved quantities: 



(4.7) 



When rp = 2 Eq. reduces to the conditions: 

0=7(-l + r)V (l + r), 
^2/''= -a + Apd^ +Qrl3d^ + 2r'^ I3d^ -^rd'^-1 

-14 d^r^-lr^d'^-i-r^d'^-i, 
0= -a + ipd^ +6f3d^r~8d^-fr + 2f3d^r^ 

-^Ud^jr^-ld'^-fr^-d'^-fr^. (4.4) 



M = 5M,P^'-^;E. 
V 11 88 

Thus again we find that the relationship: 

2c c 



25V5c27r 
176 
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(4.8) 



is exact. 



5 



To understand whether these solutions are stable solu- 
tions we use the variational principle using a trial wave 
function of compact form: 



u = Acos[dy]; < dy < ^; 



y = x- q{t) (4.9) 



where the relationship 



AdP 



pertains. (See the appendix for details) If we insert this 
into the action, we obtain for the reduced Hamiltonian 
(where L = —Pq — H[P]), using the above values of a,/3 
and 7: 



H = — (-9d^ + - 3d) 

TT 

dH 



(4.10) 



From Lagrange's equations ^ = — ^ we obtain a 
IE IP = —q = —c. Minimizing the energy with respect 
to d gives back the exact result that = 1/5 which 
is a minimum of the energy and d^ = 1/3 which is a 
maximum of the energy. This suggests that the narrower 
compacton with rf^ = 1/3 is unstable. We have confirmed 
this numerically. We have also shown numerically (see 
below) that if we start with initial compact data which 
is wider than the compacton with = 1/5 it breaks up 
into a number of these compactons. 



2. rp=4 case 

Let us now consider the case ^ = with ci = 0. 
When 7 = one of the consistency conditions Ac = 
can only be satisfied for static solutions. When 7^0 
the consistency conditions leads to : 



d^ 

AP = 



/3 



7(r^ 



6r 
c 



2)' 



(r- 1) r2 (5 + r) d'^j 
with the parameters also obeying the constraint: 
(2 + r) (4 + r) {-A + Ar + r^) (3'^ 



a = 



(-2 + 6r + r2) 7 



(4.11) 



(4.12) 



Let us look at 2 examples from this class of solutions: 

For the case 

(i) p = TO = n = 1; ^ = 

the solution is of the form: 

u = Acos'^{dy) 

Choosing for example 7 = 1/38 and /3 = 1, the above 
relations yields d = 1, A = 19c/216 and a = 672/19. 
Thus the solution is 



19c 4 
■u = — cos {x- ct) 



(4.13) 



For this solution we have that the global quantities are: 



M = 



19c7r 
"576"' 



P = 



12635c^_ 
11943936' 



E = - 



12635c^ 
17915904 



So that E = — |cP. We show in the appendix that this 
solution can also be found by minimizing the effective 
action using the compacton ansatz. In that case one also 
gets another solution with d'^ = j^, which is the maxi- 
mum of the Hamiltonian, with d^ = I being a minimum 
as a function of d. The Hamiltonian as a function of the 
width d for the above choice of parameters is given by: 



H = 



192 Vd (-77 + 19 d2 - 2 d'^) Pi 
95\/357r 



For the case 
(ii) p = m = n = 2,l~Q 
In this case we obtain the conditions: 

'^~^' " 147' ^ - 28rf4^ 



(4.14) 



Choosing /3 = 1, 7 = ^, a = we obtain c = 2A'^ and 
the solution is 



U = \ l ^ cos^ (x ■ 



ct). 



(4.15) 



For this choice of parameters we obtain for the conserved 
quantities: 



M = 



2V2' 



P = 



Sire 
"32" 



:E = — 



Sttc^ 
64 



We notice that we obtain the relationship: E/P — 
—2c/ {p+2) which we obtained by the variational method. 
For these parameters, if we use the compacton form as 
the trial wave function, and use the relation: 



A' = 



16dP 



then the effective Hamiltonian for the variational param- 
eter d is 



H 



8P' 



{lOd - 5d^ + d^) 



(4.16) 



This Hamiltonian has 2 stationary points as a function of 
d,d = 1 { minimum) and d = \f2 (maximum) The second 
solution is not a solution to the equation of motion. 

Next, let us consider a particular special solution for 
the case ci ^ 0- p = m = n = 1,Z = 0. 
For this case we have r = 2. Assuming a solution of the 
form 

u(y) = Acos^(rf?/) 
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we get the consistency equations: 



(4.17) 



which has only a trivial solution. 

For p = 2 the situation is simpler and one obtains when 
r — 1 the two relations: 



Y^[/3±(/3^-«7) 



l/2i 



A = 



8d2(/3-87d2)- 



(4.18) 



(4.19) 



For the special choice of a = 5, /3 = 3 and 7=1, there 
are two real solutions for c?^ , corresponding to (P = 1/12 
and (P = 5/12. The first solution which is stable and 
which we will discuss further in our section on numerical 
simulations is : 



9c 2. 
TT COS ( , 

14 ^ V12 



(4.20) 



This solution has for its constants of motion: 



M 



V37T 

14 

Thus we obtain: 



2437rcV3 
1568 ' 



2349 TT 
21952 ■ 



E 
P 



29c 
12" 



This shows a failure of the relationship 3.12 



The second solution which we found to be numerically 
unstable is 



u = cosily ^{x-ct)]; OO. 
The conserved quantities E, M , P are given by 



(4.21) 



97^c^/3 „ 2437^c2^/3 ^ SIttc^VS 
M = P = ^\ E ^ 
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Now we find that E/P = -c/6 ^ -2c/(p + 2). 

The special solutions with ci 7^ are not obtainable 
from a variational calculation but these are a very re- 
stricted class of solutions. 



B. Case p = m = l,n = 

For this case inserting a trial solution of the form 
u = A cos*" dy usually leads to an overdetermined set of 
equations. 

For example for p = 1, r = 2 we obtain the conditions 
(x = cos dy): 

= -ci - 48A^d^-fsm{dy)x + {Ac ~ 8A^l3d^) 

1 



2/3d^ + 5d^-f 



and 



= 18d^7 + Ud^P - a 

These equations can have 2 or even 3 positive solutions 
for d^. One particular case is a = 216, /? = 21, 7 = —2. 
In that case we get two positive solutions for d^ namely 
d = 1 leading to : (i) 



/c72 



cos(x — ct) 



(4.24) 



where the constants of motion are 
M 

and d'^ — 2 giving the solution: (ii) 



1^2;P=^;E = ~^ 
2^ ' ' 128' 256 



u — — cos \/2{x — ct) 
and we obtain for the conserved quantities: 



(4.25) 



M = Vc/S; P 



16\/2' 



E 



32\/2 



Thus both these solutions again obey the relation: 
E/P = -2c/(p + 2). The second solution with d^ ^ 2 
turns out to be numerically unstable. 
For p — 2 and r = 2 we get relations: 

ci= SA^d^j 
c- - 208A^d^-f 
-a + 48/?d2 + 1152dS = (4.26) 

as well as one constraints among the parameters 

(3 = -48d^7 

Eliminating the constraint, we obtain for the width: 



d^ = 



11527 



(4.27) 



96A^d^-fsin{dy)x'^ 



-A^a + l2A^I3d^ X* (4.23) 



As an example if we choose 7 = —3 and a — 3456 then 
we have d — 1 and for our solution: 

M= i(-^-)i/2 cos^ix- ct) 
4^39^ ^ ' 

We have not exhausted all possible solutions for this 
case, but the method for finding them should be clear to 
the reader by now. 
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C. Some Other More General Cases 

When I + n = p = m, we instead have one parameter 
family of solutions depending on the velocity c. That 
is for fixed a (3 and 7 there is a solution of different 
amplitude for different velocities c. In some special cases 
there is the possibility for two different solutions with the 
same value of c. However, in general for a given a (3 and 
7 there is only one solution with a fixed velocity c. To 
illustrate this fact, let us consider the case 

p — n; m — I — 

Inserting a trial solution of the form 

u = Acos^/Pidix - ct)]; p^2 (4.28) 

into Equation (p^ we obtain for example if 7 = l,/3 = 1, 
and a = —1/2 the conditions: 

d=|[(p+l)(2p2+4p + 3)]-i/4 



[(p+l)(2p2 + 4p + 3)]i/2 



(4-p2) 

2/3(p2 + 8p + 4) 



(4 - p^)y/{p+l){2p^ +ip + 3) 



(4.29) 



Choosing for example p — n — I, 1 = 111 = 0, "/ = 1 and 
a = —1/2 we get the single solution: 



= 8%/2cos2[ 



(4.30) 



2(18)1/4 ' 

The constants of motion for this case are: 
M = 8^(2)3/4^3; P = 4871(2)1/4^3; E = 

3v 3 

So that E/P = — ^ For the special case when p = n = 
2, one obtains 

u = A cos dy 
When TO = Z = we find 



a 
'307' 



I5O7, 



a/? 



107 



Choosing further (3 = a = \\ ^ = —1/30, we obtain 
u = VlO — c cos y 

So again for this special case we get a continuous family 
of solutions as long as c < 10. The constants of motion 
are now: 



M = 2^m^\P = -{lQ~c)-E= — (10 - c)(30 - c). 

For any particular value of p, to, n, I one can always 
find the consistency condition on the parameters a, /?, 
7, A, d, q and c so that solitary wave solutions of the 
form Acos'^[dy] exist. However we do not have a simple 
expression for these parameters for all p, to, n, I, and 
have instead in the above looked at some simple cases. 



V. NUMERICAL STUDY OF THE 
GENERALIZED KDV EQUATION 

A. Numerical Method 

In our calculations, we approximated the spatial 
derivatives with a pseudo-spectral method using the dis- 
crete Fourier transform (DFT). The equations were in- 
tegrated in time with a variable order, variable timestep 
Adams-Bashford- Moulton method. The numerical er- 
rors were monitored by varying the number of discrete 
Fourier modes between 128 and 512 and varying the es- 
timated time error per unit step between 10~^ and 10~^ 
to insure that the solutions were well converged. Mass 
and momentum were conserved to an accuracy of at least 
10~^ and the Hamiltonian was conserved to an accuracy 
of better than 10^^ in all of the calculations. 

The numerical approximation must respect the deli- 
cate balance between the nonlinear numerical dispersion 
terms in the equation. For example, when the third 
term in (2.1) is expanded, it has a diffusion-like term 
2/3m{u'^~^Ux)uxx- On the trailing edge of the compacton 
Ux > and this term acts like a destabilizing back- 
ward diffusion operator. The solution would be unsta- 
ble if it were not for the stabilizing nonlinear dispersion. 
This balance is easily lost in numerical approximation 
if the aliasing, due to the nonlinearities, is not handled 
carefully. To identify numerical artifacts due to aliasing 
and other discrete affects, we solved the equations with 
the nonlinear terms expanded in different formulations. 
For example, we compared the solutions of the equations 
when they were differenced in conservation form and non 
conservation form to identify possible numerical inaccu- 
racies. 

Also, the lack of smoothness at the edge of the com- 
pacton reduces the spectral method to first order near 
the edge and introduces dispersive errors into the calcu- 
lation. To reduce these errors and the errors due to alias- 
ing, we filtered the time derivatives by explicitly adding 
an artificial dissipation term to the equations. This term 
was defined in Fourier space to approximate the effects 
of linear second order dissipation (AxUxx) on the top 1/3 
of the Fourier modes, have no effect on the lower 1/3 of 
the modes and used a linear transition between the two 
regions. 

Most of the calculations shown here solved the conser- 
vative form of the equations with 128 DFT modes and a 
time error of lO^'^ per unit time. 



B. Numerical Investigations 

For the original Rosenau and Hyman compactons, nu- 
merical investigations 
showed some remarkable properties- namely that what- 
ever initial compact data was given, it eventually evolved 
into compactons. When two compactons scattered any 
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energy not in the original pair of compactons emerged as 
compacton-anticompacton pairs. We will find that the 
compactons of this fifth order generalized KdV equation 
have similar properties to those previously found in the 
studies of Rosenau and Hyman on their third order gen- 
eralized KdV equation. 

The first generic feature of these equations is that 
arbitrary initial compact data, as long as the width 
of the packet is larger than that of the compacton 
evolves into several compactons with the number de- 
pending on the initial energy. We show this for two 
different cases. The first case is related to the com- 
pacton of Eq. (4.20). We start of with an ini- 
tial pulse which is four times the width of the com- 
pacton and watch it evolve. This is shown in Fig. 1. 




X 



Fig. 1. Pulse with an initial width four times that of the 
compacton of Eq. (4.20) pertaining to the parameters 
p = m = n = 1,1 = and q: = 5, /3 = 3,7 = 1 namely: 
Wo = fs'^^^^ij^^)- The initial wide pulse breaks into 
compactons that collide elastically. Note the phase shift 
of the slower pulse after colliding with a faster, higher 
compacton. 

In figures 2 and 3 we show the same phenomena for the 
compacton system described by Eq. (4.7), again start- 
ing from initial data wide compared to the compacton 
solution. 




Fig. 2. Similar situation as in Fig. 1 for the parameters 
p = m = n = 2,1 = and a = 6,/3 = 4,7 = 3 relevant 
to Eq. (4.7). An initial compact wave (solid line) uq = 

cos( ^"^^° ) wider than a compacton with breaks into 

a string of compactons with the shape j4cos(^^). 




Fig. 3. Same compact wave break up as in Fig. 2 dis- 
played differently. 

In figures 4, and 5 we show similar features of the 
breakup of a compact wave for the compacton Eq. (4.13). 
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Fig. 4 Break up of a compact wave with four times the 
width of the compacton: p = m = n = 1,1 = and 
a — ^,/3 — 1,7 = An initial compact wave (sohd 
hue) uq — 2jyq cos^{ ^~J'^ ) breaks into a string of com- 
pactons with the shape ^cos''(a; — ct) by time t = 10 
(dashed hue). 




Fig. 5 Gray scale contour plot of the evolution of the 
compactons in Figs. 4. 

Then next generic feature is what happens when two 
compactons of different speeds collide. The compactons 
remain coherent and experience a phase shift. This is 
shown in Fig. 6 for the compactons described by Eq. 
(4.20). 

Unlike the solitons of the original integrable KdV equa- 
tion, these compactons after scattering also can cause 
"pair" production of a compacton and anti compacton. 
This is shown in Fig. 7 in the resolving of the wake left 
behind during the collision shown in Fig. 6. 



CO 
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Fig. 6. Collision of two solitons for the case p = m = 
n = 1,1 — and a = 5, = 3, 7 = 1. Two compactons 
described by Eq. (4.20) collide. One has speed c = 2; 
and the other has speed c = 1; Note the phase shift in the 
slower compacton after the collision. These compactons 
remained coherent, even after dozens more collisions. 



3 
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Fig. 7. A ripple (solid line) is created when the com- 
pactons first collide in Fig. 6. This ripple was extracted 
from the solution at t = 25 and used as an initial con- 
dition. By time t = 500, the ripple (dashed line) has 
separated into compactons traveling in opposite direc- 
tions. These compactons have a shape proportional to 
Acos^(^). 

When a compacton and an anti-compacton collide, as 
well as when one starts with initial data which is narrower 
than the width of the stable compacton one finds numer- 
ically blowup at late times. Whether this is an numerical 
artifact is not absolutely certain. This effect is shown in 
Fig. 8 for the same compactons as in Figs. 6 and 7. 
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T = -Pq-H[A{P),hl 



(5.2) 




Fig. 8. Possible blowup after a compacton anti- 
compacton collision, p = m = n = \,l = Q and 
a = 5,/3 = 3,7=l. Two conipactons described by 
Eq. (4.20) with speed c +1 and c = -1 collide. The 
numerical solution breaks down slightly after t — 15. It is 
not clear if the break down is due to the steep gradients 
in the solution, or because there is a true singularity that 
develops in the equations. 



APPENDIX- EXACT VARIATIONAL ANSATZ 

In this section we determine to what extent we could 
recover from the variational ansatz the exact solitary 
wave solutions we have discovered earlier by trial and 
error. That is if we assumed solutions of the form 

A{t) cos'-[6(i)(a; - qm- ^<b{x~ q{t)) < I 

for compact solitary waves and 

Ait)sech''[b{t){x- q{t))] 

for ordinary solitary waves would we recover all the exact 
solutions? 

First let us show that for the KdV equation and for the 
generalized KdV equation that we investigated earlier, we 
indeed obtain the exact solution. The Lagrangian for the 
KdV equation is 

f 1 1 

The conserved Hamiltonian is given by 



H 



dx [{ifx)^ + ^i'Pxx)'^] 



Assuming the trial wave function: 

ip.^ = u{x,t) = A{t)sech^[b{t)[x - q{t)] (5.1) 
we find that the reduced action is 



where 



and 



, 1 n , 2A^{t) 



We can rewrite H in terms of A as follows: 



5 

where we have used the two possible solutions: 



(5.3) 



(5.4) 



(5.5) 



^ 2 ' 

Since the Hamiltonian is independent of q, we have that 
P is conserved. 6 is a variable of constraint and is elimi- 
nated by the equation: 



(5.6) 



Only the negative choice for A in terms of P yields a 
positive solution for b namely: 



(6P) 



1/3 



42/3 



(5.7) 



Eliminating b the reduced action is 



(5.8) 



Varying the action we find the velocity is a constant: 
q = (3P)2/3 = c 

Thus A = — I and we get the usual exact answer: 



u{x, t) — — -sech [—^{x — ct)] 



We also find that 



H = -P{-{QbPY'^ + b^) 
5 



(5.9) 



(5.10) 



has a minimum at the exact value of b for fixed P. 

Next we consider the class of exact compact solitary 
waves that we found for the generalized KdV equation of 
ref. P] . In this case the Lagrangian is 



f 1 1 

L = J dx i-^^^t + a-^—^{^. 



and the Hamiltonian is: 
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Now we assume a solution of the form (r = 2/p) 

Lp^ ^u{x,t) ^ Acos^/P[d{t){x - q{t))] (5.11) 
We obtain for the reduced action 

T = -Pq-H[A,d] (5.12) 

where now 

A2v^r(i/2 + 2/p) 



p = 



2d{t)T{l + 2/p) 



(5.13) 



H = 



A^+P{AI3d\2 + 3p + p2) _ a(4p + p2)}0Fr(i + |) 



2dp2(i+p) (2+p)r(2+|) 



(5.14) 



Replacing A by 



/ 2dFr(l + 2/p) \ 

\^ni/2+2/p)) 



1/2 



we obtain: 



i? = (2rf)tpi+f 

{4/3^2(2 + 3p + p2) _ a(4p + p^)} rf (1 + f ) 



£t^£/1 , 2n (5.15) 



(2p + 5p2 + 4393+p4) 7,4r2(i + ^) 



We determine the constraint variable d by ^ — , 
we obtain 



ap 



4/3(p+l)(p + 2) 
Lagrange's equations gives: 

dH P+2H 



(5.16) 



(5.17) 



We then have that 



AP 



c(p+l)(p + 2) 
2a 



which has a minimum at cZ^ = 1/12. 

Now let us look at our generalized equation when 7 7^ 
0. For the special case p — m — n;l — considered in 
this paper we have that the Lagrangian is: 



L{p ^ m ^ n;l ^ 0) ^ [ dx [^ipx^Pt + a- — [^'^ 
J 2 (p + 1 



p+2 



ip+l){p + 2) 

- Piv.n^..)^ + lylvL^], (5.19) 

Introducing a trial variational function of the form: 

u = Aeon'' d{t){x - q{t)) (5.20) 
with r ^ A/p and the constraint: 

_ (2 + r)(4 + r)(-4 + 4r + r^)/?^ 



7(-2 + r2 + 6r)2 
Using the fact that 



(5.21) 



P - A^^^ 



2 ^r(l/2 + r) 



2d r(l + r) 



to eliminate A in favor of P, we again find we can write 
the reduced action as 

j dt{^Pq- H[P,d]} (5.22) 

where 

H= (2rf)2/'-pi+2/v2r[l + r]i+2/'- X 

(12/32 + spd^j - 32d^j^ + 20P^r - S2f3d'^jr + 176d^-/'^r 
-19/3^2 - 32(3d^jr^ - 172dS^r^ - 24^^^ 
+152pd'^^r^ - 152d'*7V3 - 4/3^"* + 50pd^-/r'^ 
-152dSV'' + APd^jr^ - 40d''7V^ - 3d'^-f^r^) x 

(4^^i/r^_2 + 6r + r^)\[3 + r]T[l/2 + rf^'')-^ (5.23) 

From the equation that eliminates the constraint vari- 
able d: 



dH 
~dd 







and recover our previous exact result [^: 

(5.18) 

As a function of d for fixed P, _ff is a minimum at the 
constraint equation value of d. As an example, when 
p = 1,P = 1 and (3 ~ 1/2, a = 1 one obtains for H[d]: 

H^l{^)'/'i-5 + 12d') 



we find there are two solutions for d^. One solution 



d' = 



7(7-2 + 6r- 2) 



(5.24) 



is a minimum of H [d] for fixed P and is an exact solution 
of the generalized KdV equation. The other solution for 

d2 



13 (3 + 2r) (-4 + 4r + r2) 
7 (-2 + 6r + r2) (8 + 4r + 3r2) 



(5.25) 



is a maximum of the energy H[d] for fixed P and is not a 
solution of the equation of motion. An example discussed 
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earlier is the case p = r = 2 with /? = 1, 7 = 1/14, and 
q; = 96/7. In that case we have 



op2 

H[d] = ——{10d-5d^ + d^), (5.26) 

with two extrema: d = 1 which is a minimum of H 
and yields an exact solution, u = y/c/2cos^{x — ct), 
and d = \/2 which is a maximum and leads to u = 
\/3c/2 cos^ V2{x — ct) which is not a solution of the orig- 
inal generalized KdV equation. 
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